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In this supplementary material, we investigate further the impurity-induced freezing mechanism 
in a doped system of 3D weakly coupled ladders resembling BitCui^Zn^ZnPOg using QMC. 



I. INTRODUCTION 



In the previous letter [1], we have argued that the collective freezing of effective moments having a 3D 
extension V£ ~ £,x£,y£,z at T > T g is actually controlled by the exponentially decaying 3D coupling of the 
general form 

|Jfg(^ J 3D exp(-|-|--£), (1) 

expected to occur for the wide class of spin gapped materials [2, 3, 4, 5, 6, 7]. The average coupling J avg 
taken over all possible J^%{r) does account for the broad distribution of effective interactions and is just 
given by 

</av g = <|J 3 C ^1l>^3D— ^ (2) 

where ~ £,x£,yiz is the magnetic volume occupied by each induced moment. We propose that this average 
coupling governs the ordering, i.e. T g ~ J avg . 



II. MICROSCOPIC MODEL 



We want to check such an analysis against QMC simulations on a diluted 3D model of weakly coupled 
ladders (schematized in Fig. 1) with the following parameters: J±/ J = 0.1 and J^d/ J = 0.05 which, 
using a value of J ~ 100 K corresponds to a spin gap A ~ 35 K and a transverse 3D coupling J^o ~ 5 
K. We then introduced non-magnetic impurities (open circles in Fig. 1) and performed large scale QMC 
simulations on 3D samples of sizes L x L x L/2, with L = 16, 24, 32, 48, down to temperature T / J = 0.01. 
We also performed disorder averaging over a large number N of independent disordered samples ranging 
from N = 500 for L = 16 to N = 100 for L = 48. 
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FIG. 1: (Color online) Schematic picture of the 3D system of coupled ladders used for the QMC simulations. 

III. QMC RESULTS FOR THE CRITICAL TEMPERATURE 

The results for the 3D ordering temperature T g are shown in Fig. 2 versus the impurity concentration x. 
The transition was found by the standard technique using the finite size scaling of the spin stiffness at the 
transition point, as examplified in Fig. 3 and discussed below. As shown in Fig. 2, we get a linear increase 
for T 9 (x) up to a threshold ~ 3% where T g saturates as expected. The linear part can be fitted by the form 
T g = 2.8x which compares quite well to our estimate Eq. (2). Indeed, with our parameters, we expect an 
average coupling J avg pa J^n x x x pa 0.05 x 30 x x = 1.5x, meaning that with such a definition for 
V£, we get T 9 (x) w 2 J avg . 



IV. DETAILS ABOUT THE CRITICAL POINT 



The way the critical ordering temperature was extracted from QMC simulations on finite size systems 
is actually standard since it relies on the finite size scaling of the order parameters, as for instance used 
in Ref. [8]. Therefore we computed the spin stiffness p s , directly related to the square of the AF order 
parameter. In an 3D AF ordered phase p s is finite whereas it is in a disordered phase. At the critical point 
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FIG. 2: (Color online) Critical temperature TG (in units of J) plotted versus impurity concentration x. 
between the two regimes, there is a well-known finite size scaling 

Ps (L) ~ L 2 - D - z , (3) 

where D is the space dimension (here D = 3) and z is the dynamical exponent (z = for a finite temper- 
ature phase transition). Therefore we expect p s x L to be a constant at the critical point where a crossing 
of the various system sizes should occur. We thus used such a criterion to identify the ordering transition 
atx = 1, 2, 3, 4, 5%. Results of such an analysis for a concentration x = 2% are displayed in Fig. 3 
for L = 16,24,32,48 with a critical point found at T g = 0.055 J. In Fig. 3(A), we show the average 
stiffness versus Tj J. In fact the spin stiffness is a directionnal quantity and can thus be computed in all 
space directions x, y, z, or averaged over all directions. The crossing of p s x L is shown in Fig. 3 (B) as 
well as in insets (X,Y,Z) for all the components of the stiffness. This clearly shows that the ordering is 
fully three dimensional and we find a remarkable agreement for the crossing temperatures in all directions 
at Tg/J = 0.055. 
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FIG. 3: (Color online) QMC results for the 3D model of coupled ladders with x = 2% of non-magnetic impurities. 
(A) shows the average 3D stiffness versus Tj J. The same data are shown in (B) where p s x L has a crossing for all 
sizes at T/J = 0.055. Insets (X,Y,Z) show all the components of the spin stiffness times L that cross at the same 
critical temperature. 

V. COMPARISON TO EXPERIMENTS AND CONCLUSIONS 

To finally conclude on this issue of the 3D transition, we carefully checked that the ordering transition 
is a true AF 3D ordering which occurs at a freezing temperature T g proportionnal to the average coupling 
as proposed in the paper [1]. We indeed confirm a linear regime with x at low concentration followed by a 
saturation a larger x corresponding to the fact that the average distance between impurity start to be of the 
order of the correlation length £. As a comparison we plotted on a common graph (Fig. 4) experimental 
results for T g rescaled to their x = 3% values for various spin-gapped materials together with the QMC 
results of this study. The agreement is very good. 
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FIG. 4: (Color online) Transition temperatures T* (rescaled to their values at x=3%) versus impurity concentration 
for various low-D spin-gapped systems: coupled ladders Bi(Cui_ x (Zn or Ni^^POe from this study; isolated ladder 
Sr(Cui_ a; (Zn or Ni^^Os; Haldane chain PMNii-zMg^'V^Os; spin-Peierls chains Cui-^Zn or Ni^GeOa. QMC 
data of Fig. 2 are also shown for comparison. 
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